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We show that a combination of a half-cycle pulse and a short nonresonant laser pulse produces a 
strongly enhanced postpulse orientation. Robust transients that display both efficient and long-lived 
orientation are obtained. The mechanism is analyzed in terms of optimal oriented target states in 
finite Hilbert subspaces and shows that hybrid pulses can prove useful for other control issues. 
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Laser controlled processes such as molecular alignment 
and orientation are challenging issues that have received 
considerable attention both theoretically and experimen- 
tally |l|. Whereas strong nonresonant adiabatic pulses 
can exhibit efficient alignment and orientation only when 
the pulse is on 0, 3 , linear polar molecules can be ori- 
ented under field-free conditions after the extinction of a 
short half-cycle pulse (HCP) pi, @. Its highly asymmet- 
rical temporal shape imparts a sudden momentum kick 
through the permanent dipole moment of the molecule, 
which orients it. This extends the use of a perma- 
nent static field (combined with pulsed nonresonant laser 
fields) 0,13- Similar lyto the alignment process by a non- 
resonant short pulse |3| (measured by (cos 2 0) with 9 the 
angle between the axis of the molecule and the polarisa- 
tion direction of the laser field), the orientation by a HCP 
(measured by {cos 9)) increases as a function of the field 
amplitude until it reaches a saturation at (cos#) « 0.75, 
which corresponds to an angle cos -1 (cos 0) « 41°. Over- 
coming this saturation has been theoretically proved with 
the use of trains of laser pulses (HCP kicks) in the case 
of alignment |9j (orientation |l0j). For applications it 
is of importance to reach an efficient orientation. An- 
other crucial point, that has received less attention so far 
|ToL ITl| , is the duration during which the orientation is 
above a given threshold, which one would like to keep 
as long as possible. An important step made in was 
to establish a priori the two oriented target states (of 
opposite direction) in a given finite subspace generated 
by the lowest rotational states. These target states are 
optimal in the sense that they lead respectively to the 
maximum and minimum values of (cos 9) in this given 
subspace. The choice of a suitable small dimension of 
the subspace allows one to generate an oriented target 
state of relatively long duration. The identification of 
such an optimal target state opens in particular the pos- 
sibility to use standard optimization procedures (see e.g. 

E2)- 
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One of the main challenges consists now in reaching 
such an optimal target state in a subspace of low dimen- 
sion, characterizing a long-lived and efficient orientation, 
by a simple external field in a robust way and to ensure 
the persistence of this effect with respect to thermal av- 
eraging of finite temperature. We propose in this Letter 
a process that possesses such properties. By superim- 
posing a pump laser field to a half-cycle pulse, we show 
that the maximal orientation reached after the pulse is 
significantly beyond the one induced by an HCP, and 
that it displays furthermore a larger duration. We ob- 
tain in particular the saturation (cos#) « 0.89, which 
corresponds to the angle cos -1 (cos 9) ~ 27°. This effi- 
cient and long-lived orientation is obtained by adjusting 
only two parameters: the amplitudes of the laser and 
HCP fields. We obtain robust regions of the parameters 
generating this orientation. 

We show that this process allows one to approach an 
optimal target state in one step. Under the action of 
a single hybrid pulse, the number of rotational states 
that are significantly populated remains finite and con- 
trollable. This generates a finite dimensional subspace in 
which an optimal target state can be constructed. When 
the dimension of this subspace increases, the associated 
optimal state yields a higher orientation efficiency while 
its duration decreases. By choosing appropriate intensi- 
ties of the pump laser field and of the half-cycle pulse, 
we can both select and reach the target state with the 
desired efficiency and duration. 

We consider a linear molecule in its ground vibronic 
state described in the 3D rigid rotor approximation. The 
effective Hamiltonian including its interaction with a 
HCP simultaneously combined with a pump laser field 
of respective amplitudes £hcp(£) and £h(t) is given by 

H cS (t) = BJ 2 - o H op(*) cos 9 - <z L (i) cos 2 9, (1) 

where B is the rotational constant, ajjcp = Mo £hcp with 
Ho the permanent dipole moment, and a-^ — Aa£ 2 /4 
with Aa the polarisability anisotropy. Note that Acv 
is positive for linear molecules, which gives positive val- 
ues for a,L whereas the sign of <zhcp is determined by 
the sign of the HCP amplitude £hcp- The dynamics of 
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the system is readily determined with the help of the 
propagator in the impulsive regime, where the duration 
t of the pulse is much smaller than the rotational pe- 
riod T rot = nh/B For the process we suggest here, 
in the dimensionless time s = i/r ro t whose origin coin- 
cides with the extinction time of the pulse, the propaga- 
tor reads U(s,0) = e - Mj2s e lAHCP cos V Al cos2 e . The pa- 
rameters Ajjcp = \ J dt cihcp {t) and = j- J dt (t) 
are respectively the total dimensionless areas of the HCP 
amplitude and of the laser intensity. 

We first consider the case of a cold molecule whose 
state after the pulse reads \<f>(s)) — U(s,Q)\j = 0} where 
\j) stands for the spherical harmonics Y™ with m = 0. 
The orientation is measured by the expectation value 
(cos0)(s) = (<j)(s)\ cos 8\<f>(s)) . It is well-known that it 
is a periodic function (of unit period), which exhibits 
peaked revivals corresponding to molecular orientation 
along the field direction 0, ||J . 
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FIG. 1: (Color online) Numerical contour plot of 
maxjj | (cos 6) | as a function of the total (dimensionless) areas 
Al and Ahcp- The straight line Ahcp = 2.5 Al approxi- 
mately indicates the maximum for given Ahcp- 

The maximal field-free orientation is displayed in Fig. 
H by plotting the maximum of | (cos 9) | reached over a 
period as a function of the total areas A L and A H cp at 
zero rotational temperature, calculated with the above 
propagator. The case of a single HCP coincides with 
the ordinate axis. In the absence of the laser pulse it is 
seen that the maximal orientation saturates to a value of 
(cos0) w 0.75. In the presence of a simultaneous pump 
laser pulse, one observes a wide two-dimensional plateau 
as well as an island which are both associated with an 
orientation much higher than the saturation limit of a 
single high intensity HCP. The plateau is relatively flat 
in a large two-dimensional region centered approximately 
around the line ^4hcp = 2.5 A^, implying that one can 
robustly reach a high efficiency for the orientation with 
a moderate HCP intensity. The island observed around 



^4hcp = 1-25 and Al = 3.7 indicates that it is also pos- 
sible to overcome the above saturation by combining an 
HCP of intensity slightly above unity with a laser pulse 
of high intensity. 

The direction of the orientation can be chosen by the 
sign of the amplitude of the HCP. Expressing the observ- 
able in terms of the projections Cj = (j\<j>(0)) of the wave 
function right after the pulse onto the rotational states 
\j) leads to (cos0) (s) = § £ c*c j+1 e- 2 ™^ +1)s +c.c. The 
coefficients Cj can be calculated for the above propaga- 
tor in the approximation (j\ cos 9\j ±1) ~ 1/2 which is 
more accurate for j > (one has (O|cos0|l) ~ 0.58, 
(l|cos0|2) « 0.52, (2|cos0|3) « 0.51, •••)■ This shows 
that the sign of each product c*Cj+i , and hence of (cos 9) , 
changes with the sign of j4hcp- For positive AhcPi the 
expectation value corresponding to max s | (cos 9) \ is posi- 
tive on the island and negative on the plateau. We show 
below that the high orientation efficiency obtained in the 
region along the line Ahcp = 2.5 A^ has the remarkable 
property to approach very closely an optimal state as de- 
fined in [T(ij | , which combines orientation of high efficiency 
and of long duration. 
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HCP 

FIG. 2: As a function of ^4hcp, for Al — ^4hcp/2.5, upper 
panel: (i) max s |(cos#}| (solid line, right axis);(h) for compari- 
son, same quantity with Al ~ (dotted line); (iii) duration A 
for which |{cos#)j > 0.5 (dashed line, left axis); Lower panel: 

square modulus Pn = \{X- l^(smax)) | 2 of the projection of 
the state at the time giving the maximum of |(cos#)| on the 
optimal state |x_ ) (solid line, right axis); dimension iV of 
the corresponding subspace (step function, left axis). 

The upper frame of Fig. [^compares the maximal value 
of | (cos 0)| as a function of ^hcp (i) without the laser, 
and (ii) along the straight line Ahcp = 2.5 A^ seen in 
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Fig. 1. This shows that the saturation of (cos 9) w 0.75 
obtained with the HCP alone is significantly overcome 
(up to (cos 9) w 0.89) when the HCP is associated with 
the laser of appropriate area. The upper frame of Fig. [21 
also illustrates the duration of the orientation defined as 
the time during which |(cos0)| > 0.5 (see also Fig. OjJ. 
It is seen that the duration of the revival can be as large 
as about 18% of the rotational period. The properties 
displayed in Fig. 2 are robust with respect to the pa- 
rameters Ahcp and which need not be in a strict 2.5 
ratio. 

The efficiency of the obtained oriented state and its 
large duration are explained in terms of an optimal state 
as defined in [Tfij . We recall that the optimal states 
correspond to the two states that respectively minimize 
and maximize the projection of cos9 in the finite sub- 
space Hn spanned by the N lowest rotational states 
|0),|1) ••• ,\N-1), namely cos (Ar) 9 = U N cos 9U N with 
the projector IIjv = X^o* U)0l- Considering a finite 
subspace yields an operator that has a discrete spec- 
trum, whose eigenvectors are readily calculated and for 
which the duration of the orientation provided by these 
states can be computed. Furthermore, the controllabil- 
ity of the system can be completely analyzed 0, Il5j . 
For a given dimension N, the two optimal states are the 
eigenvectors associated respectively with the smallest and 
the largest eigenvalues of cos^' 9. In the approximation 
(j\ cos 9\j ±1) ~ 1/2 one obtains 



[lCj . The lower frame of Fig. [JJshows that the HCP-laser 
combination of appropriate areas leads in a single step 
to a wave function that is remarkably close to the opti- 
mal state Ix^) (more than 90% for 1.5 < A H cp < 5). 
Figure [5] also indicates how the dimension N of the em- 
bedding subspace can be chosen by the value of j4hcp 
with j4l = ^hcp/2-5. Notice the linear character of 
this necessarily stepwise function. In the island region of 
Fig. QJ, the dynamics also generates a state close to an 

optimal one: Kxfl^w))! 2 w 0.86 for ^ HC p = 1-25 
and Al = 3.7. For comparison, we note that the same 
optimal states are reached in [l(j by a different process 
involving 15 short HCP kicks sent at specific times and 
with a low amplitude in order to remain in a given sub- 
space. 




\X± > 



■ JV-l 



N- 



3=0 



tE^'+M^Ti)^ (2) 



giving the approximate optimal orientation 
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The (relative) duration A of the orientation is defined 
as the time during which | (cos 9) | > 7 for the revival 
of maximum efficiency, with 7 arbitrarily chosen as 1/2. 

We can determine the duration Ajv for the state |x± ) 
by summing the above expression for (cos 9) (s) and ex- 
panding the result to second order around its extremum, 
obtaining 




1 — 7/cos 



N+l 



(4) 



where T N = a(N + l) 2 - (N + l) with a = 2/3 - 1/tt 2 . 
The shape of Ajv as a function of N is similar to the 
dashed curve on the upper panel of Fig. 2, independently 
of the specific value of 7. In particular, the decrease of 
this duration for large N is due to the factor Tn- 

Finding a process that drives the system to an op- 
timal state x±^ guarantees an efficient orientation to- 
gether with a large duration if the dimension N of the 
subspace TLm generated by the dynamics is relatively low 



FIG. 3: Orientation as a function of dimensionless time for 
(i) j4hcp = 3 and A-l = 3/2.5 (solid line) and (ii) the optimal 
state X- with a time translation (dashed line). 



Figure 3 shows an example of orientation, measured by 
(cos 9) , as a function of time for a point on the straight 

line of Fig. 1 (A H cp = 3). The result is close to that 

(5) 

given by the optimal state x~ whose minimum is taken 
at the minimum of (cos 9) generated by the hybrid pulse. 

We obtain in this case (x*? |0(s ma x)}| 2 ~ 0.93 (with 
less than 2% of the total population outside the sub- 
space H5.) The maximum orientation (in absolute value) 
I (cos 6) I Ri 0.9 occurs at s ~ 0.9. One can observe a rel- 
atively large duration of the orientation. In contrast to 
the case of a sole HCP, the presence of the laser pulse of 
appropriate area (A^ = Ahcp/2-5) allows us to obtain, 
immediately after the pulse, projections cj on the rota- 
tional states whose moduli are very close to the moduli 
of the corresponding components of both optimal states. 
The phases of these projections just after the pulse gen- 
erally differ from those of the components of the optimal 
states, but are brought by the free evolution closer to 
those of one or the other optimal state. In the case of 
the plateau region, the set of phases after the pulse leads 
to the state x^ f° r positive Ahcp and yields thus a 
minimal value for (cos 9) . This analysis extends to the 
island region where revivals of opposite sign are observed. 
As discussed above, by changing the sign of Ahcp while 
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keeping fixed, one obtains orientation revivals of the 
same absolute values but opposite sign. The direction of 
the orientation can thus be controlled by the sign of the 
HCP pulse. 
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FIG. 4: Same as Fig. Q for the dimensionless temperature 
f = 5. 

Considering the effect of temperature amounts to sta- 
tistically average over the solutions of the Schrodinger 
equation with different initial conditions \j, m) weighed 
by a Boltzmann factor. Figure 0] shows the maximal ori- 
entation, measured by the appropriate expectation value 
of cos 0, as a function of the field parameters for a dimen- 
sionless temperature T = kT/B = 5 (which corresponds 
to T « 5 K for the LiCl molecule). Notice that the is- 
land disappears while the region around the straight line 
Al = Ahcp/2-5 persists. The efficient and long-lived 
orientation revivals are therefore robust with respect to 
thermal averaging and to the field parameters. The effi- 



ciency is lower than at T = OK for the same field ampli- 
tudes, but one can recover the same value by increasing 
the amplitudes along the straight line. 

In conclusion, we have shown that a combination of 
a half-cycle pulse and a short nonresonant laser pulse of 
appropriate amplitudes leads to efficient and long-lived 
revivals of orientation beyond the known saturation. Fur- 
thermore, this is achieved in a controllable manner since 
the desired target state can be chosen in a set of optimal 
target states defined in Hilbert subspaces of low dimen- 
sion and be reached with a projection larger than 90% 
by a single hybrid pulse. As an illustration, the ground 
state of a KC1 molecule with rotational constant B 0.13 
cm -1 (T ro t ~ 128 ps) and dipole moment hq « 10.3 D 
gives ^4hcp ~ 3 for a pulse duration of 2 ps and a HCP 
amplitude of 100 kV/cm. To be on the optimal line of 
Figs. [I]or 4 requires A\, w 1.2 which corresponds to a 
peak intensity / m 10 11 W/cm 2 for the laser field. These 
parameters lead to max s |(cos0)| w 0.85 an a duration of 
approximately l/8 th of the rotational period for a cold 
molecule (see Fig. [3J|, and to max s |(cos0)| ~ 0.73 and a 
duration of approximately l/20 th of the rotational period 
for T = 5 K. The interest of hybrid pulses is not limited to 
molecular orientation but extends to optimization issues 
of a large class of systems where symmetries need to be 
broken or selectively addressed (e. g. the control of tun- 
neling). The central element consists in using an external 
field that plays individually on couplings of different sym- 
metries. In order to drive the dynamics even closer to an 
optimal target state, standard optimization algorithms 
can be used for trains of these hybrid pulses (with for 
instance the delays and/or the relative amplitudes be- 
tween the kicks, or even a delay between the HCP and 
laser pulses) and should require only a low number of hy- 
brid pulses since the first step already brings the system 
very close to the target state. 
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